A two-dimensional analytical modeling for channel potential and threshold voltage of short channel triple material symmetrical gate Stack (TMGS) DG-MOSFET
Tripathi Shweta†,
Department of Electronics & Communication Engineering, Motilal Nehru National Institute of Technology, Allahabad-211004, India

 

† Corresponding author. E-mail: shtri@mnnit.ac.in

Abstract
Abstract

In the present work, a two-dimensional (2D) analytical framework of triple material symmetrical gate stack (TMGS) DG-MOSFET is presented in order to subdue the short channel effects. A lightly doped channel along with triple material gate having different work functions and symmetrical gate stack structure, showcases substantial betterment in quashing short channel effects to a good extent. The device functioning amends in terms of improved exemption to threshold voltage roll-off, thereby suppressing the short channel effects. The encroachments of respective device arguments on the threshold voltage of the proposed structure are examined in detail. The significant outcomes are compared with the numerical simulation data obtained by using 2D ATLAS™ device simulator to affirm and formalize the proposed device structure.

1. Introduction

During the last five decades since the invention of the MOSFETs, we have seen rapid and stiff advancement in the evolution of integrated circuit engineering.[13] The principal driving force behind this tremendous development has been the downscaling of the dimensions of the MOSFET.[3,4] The greater the downscale of MOSFET, the higher its packing density becomes, the faster its operation and the lower its power dissipation is.[5] The perennial downscaling slipstream of MOSFETs has steered eminent execution and low power functioning as prescribed by ITRS, however some newfangled effects in the device performance have also arisen.[46] In particular, a gathering of undesirable sources of difficulties springs up, which are jointly called short channel effects (SCEs) that hinder further advancement in transistor downscaling.[6] In order to sustain the rate of advancement in device performance with uninterrupted downscaling, alterations to device design and fabrication are the essential requirements. In connection to this, a number of new extended device structures already have been reported in the literature to handle the severe challenges associated with SCEs as well as device efficiency.[715] One of the predominant anticipated solutions in the nanoscale regime to inhibit these SCEs is DG-MOSFET.[16] However to cope with the device downscaling and quashing the SCEs, there arises the need to allow some relatively thicker gate oxide along with high-k effectuation to exert the electrostatic integrity of the DG MOSFET device.[7,15] Due to this progression, DG-MOSFET is anticipated to be further investigated to subdue the SCEs.

A multimaterial-based gate structure is anticipated to provide further improvement in the performance of DG MOSFET since it may create a step potential profile to screen the effect of the drain on the channel.[12,13] Another possible solution to suppress SCEs is by using gate stack. If aforementioned solutions are combined in order to improve the performance of DG MOSFET then it may prove to be an effective option over the existing technology. There exist several reports that establish the superiority of the performance of TM DG MOSFET and gate stack structures.[1215] Razavi and Orouji[12] proposed TM-DG MOSFET structure and through numerical simulation results they showed that SCEs such as drain-induced barrier lowering (DIBL) and hot carrier effect are reduced significantly in comparison with conventional DG MOSFET. They also demonstrated that TM-DG MOSFET leads simultaneously to the enhancement of transconductance and the reduction of drain conductance which itself leads to higher DC gain than the classical DG MOSFET. Pradhan et al.[13] reported a simulation-based study for the performance differences among six dissimilar gate stacked device structures by utilizing gate and channel engineering. They also considered the gate stack TM DG MOSFET for the simulation purpose and established that TM DG MOSFET with gate stack structures improves DIBL significantly. Alamgir et al.[14] investigated the performance of DG-MOSFET with high-k stack on both the top and bottom gates using numerical simulation results. They observed that higher drain current in the low and high drain voltage regions, as well as lower threshold voltage and improved subthreshold swing can be yielded by considering their proposed structure. Jin et al.[15] presented an analytical-model-based study for surface potential, threshold voltage, subthreshold current and subthreshold swing of symmetrical gate stack dual gate strained silicon MOSFETs. They showed that the high-k region on the oxide layer exhibits reduced SCE.

Therefore, based on the literature it can be concluded that pertaining to the scaling down of the device thickness, triple material gate stack structure is anticipated to be a more effective option over existing DG MOSFET device. To the best of my knowledge, so far in the literature analytical modeling and simulation-based study of TM DG MOSFET with gate stack have not been done. Hence in the present paper, combining the advantages of high-k gate stack and triple material structure, a new compact analytical model for triple material symmetrical gate stack (TMGS)-DG MOSFET is presented. In this proposed device structure, the front and back gate electrode each comprising three materials with different work functions and symmetrical high-k gate stack are considered for developing the compact analytical model and numerical simulation. Moreover, the channel of the simulated model is adopted to be lightly doped in order to increase the mobility of the charge carriers and stave off the problem of threshold voltage variation due to random fluctuation of dopant atoms. The result analysis is based on the solution of two dimensional (2D) Poisson’s equation with the appropriate limit specifications employing 2D ATLAS™ device simulation tool.[17]

2. Device structure

The schematic device structure of TMGS-DG MOSFET used for analysis and simulation is shown in Fig. 1, where L (L = L1 + L2 + L3), tsi, tox, and tk are gate length, channel thickness, thickness of gate oxide and thickness of material with large dielectric constant, respectively. The x and y axes of the schematic structure are considered to be along the channel length and channel thickness, respectively.

Fig. 1. Schematic diagram of TMGS DG MOSFET, where L (L = L1 + L2 + L3), tsi, tox, and tk are gate length, thickness of silicon channel, thickness of gate oxide, and thickness of high dielectric constant material, respectively.

The gate electrode of the DG MOSFET is considered to be made up of three materials having different work functions ϕm1, ϕm2, and ϕm3 (such that ϕm1 > ϕm2 > ϕm3) and deposited over the respective lengths L1, L2, and L3 on the gate oxide layer. The gate material placed at the source end with the highest work-function i.e., ϕm1 forms a control gate, the middle material with intermediate work-function value ϕm2 forms a first screen gate, and the remaining material with the lowest work function ϕm3 placed at the drain side of the gate forms a second screen gate. In the present work tungsten disilicide (WSi2, ϕm1 = 4.8 eV), Hf0.27Ta0.58N0.15 (ϕm2=4.6 eV), and Hf0.40Ta0.46N0.14 (ϕm3 =4.4 eV) are used as the materials for control gate, first and second screen gates, respectively. The gate terminal of the TMGS-DG MOSFET device is made by tying together control and screen gates in order to feed the device with the same gate voltage Vgs. In the proposed structure, high-k gate stack is used in order to combat the SCEs therefore in place of tox effective gate oxide thickness (teff) is used. teff can be given as[15]

where tk is the high dielectric material thickness and εk is dielectric constant of the high-k material.

3. Model derivation
3.1. 2D potential modeling

The 2D Poisson’s equation for potential distribution in the lightly doped channel region can be given as[15]

where n = 1, 2, and 3 correspond to the regions under metal M1, M2, and M3, respectively. ϕn (x,y) is the electrostatic potential in the channel region, q is the electronic charge on electron, Na is the acceptor concentration of channel, and εsi is the permittivity of silicon.

The Poisson’s equation can be solved using the following boundary conditions.

The potential distribution in the channel is assumed to be parabolic third order polynomial and can be expressed as[18]

where cn1 and cn2 each are an arbitrary function of x. The Poisson’s equation is solved separately under the three metal gates by using the boundary conditions from Eqs. (3)–(13), to obtain the coefficients ϕSn (x), cn1 (x), and cn2 (x) of Eq. (15) and the values obtained can be written as

Now, substituting y = tsi/2 into Eq. (15) and utilizing Eqs. (16) and (17), we can obtain ϕSn (x) as

where ϕCn(x) is the mid channel potential of the device given by

The characteristic length associated with mid channel potential of the channel is significant in determining SCEs of the DG MOSFET.[19] Therefore, characteristic length can be determined by solving Poisson’s equation in the middle of the silicon channel.[19] Now, following this argument Eq. (2) can be solved at the center of the silicon channel as,[20]

where λ is the characteristic length associated with the center channel potential given as[17]

The generalized solution of Eq. (20) can be written as[17]

where , and

The coefficients An’s and Bn’s of the generalized solution Eq. (22) corresponding to the channel regions I, II, and III can be obtained by solving Eq. (20) in conjunction with the boundary conditions of Eqs. (6)–(13) and can be given by

Now using Eqs. (15), (18), (22), and Eqs. (24)–(29), the 2D channel potentials for the three regions of the channel can be obtained as

where n = 1,2,3.

3.2. Threshold voltage modeling

Under threshold condition, the channel under the control gate will not conduct until channels under both the screen gates are turned on. The conduction from source to drain will take place only when all the three gates are turned on. This suggests that threshold voltage of the device under consideration will systematically be governed by control gate region of length L1 with highest work function ϕm1. The effect of drain-to-source voltage on the subthreshold behavior of the device can be monitored by position and value of minimum channel potential. The position of minimum channel potential xmin can be evaluated by solving equation

and given as,

Now, minimum value of channel potential ϕC1 (x) = ϕC1 (xmin) = ϕC1 min, i.e., minimum channel potential can be obtained by using Eq. (22), and equation (31) can be written as

Now, the virtual cathode potential ϕvc (y) can be given as[21]

The inversion charge density (Qinv) depends exponentially on the surface potential. Now the inversion carrier charge sheet density at ϕvc (y) can be given as[22]

In order to obtain the analytical expression of Qinv we may use the property of undoped DG MOSFET that the effective conduction path is located at y = tsi/4 from the silicon surface. Thus, Qinv can be approximated as

Now using the value of ϕvc (y) from Eq. (33), in Eq. (35), Qinv can be given as

Now assuming Qin = QTH under threshold condition Eq. (36)[22,23] becomes

where

To obtain the threshold voltage of long channel device, the limiting condition of L → ∞ in Eq. (37) yields

Now using Eq. (37), equation (38) can be solved to obtain short channel threshold voltage of the device under consideration as follows:

4. Results and discussion

In this section, we present some theoretical results obtained by the proposed model that is validated by comparing the theoretical results with the numerical simulation data obtained by using the ATLAS™ device simulator. The drift-diffusion model has been used for carrier transport, and classical Fermi–Dirac statistics model is used for the numerical simulation.[24] To integrate the effect of high electric field in the device channel, the field dependent mobility model has been implemented in the ATLAS™ simulation.[25,26] Further, the ratio of lengths of regions I, II, and III is taken to be L1 : L2 : L3: 1 : 1 : 1 for all the computations. The p-type substrate concentration is taken to be uniformly doped with a doping concentration of Na = 1 × 1016 cm−3 in all the three regions under different gate metal contacts. The source and drain regions are both assumed to be heavily doped with doping concentration of NS/D = 1 × 1020 cm−3. The dielectric constant value of high-k material is taken to be 20. The typical values of tsi, tox, and tk are taken to be 20 nm, 2 nm, and 4 nm respectively.

Figure 2 shows the variations of surface potential of TMGS-DG MOSFET along the gate length for different values of drain-to-source voltage (Vds). It can be clearly seen from the figure that the area under the control gate is effectively screened from the variations in Vds, thereby giving rise to negligible DIBL effect in the device due to the presence of the two screening gates. Further, it may also be observed from the figure that the position of minimum channel potential lies under the control gate, which justifies the assumption that we have made for modeling the threshold voltage. The variations of threshold voltage with gate length for different values of Vds are plotted in Fig. 3. It can be seen from the figure that for longer gate length devices, threshold voltage varies little with the drain bias compared with the corresponding shift of threshold voltage for shorter gate length devices.

Fig. 2. Variations of channel potential with position along the channel length for different values of drain-to-source voltage (Vds).
Fig. 3. Variations of threshold voltage with gate length (L) for different values of drain-to-source voltage (Vds).

Figure 4 demonstrates the variations of threshold voltage with gate length for different values of silicon thickness (tsi). It is noteworthy here that for the lower values of tsi, the roll-off in the threshold voltage is relatively low in comparison with that for the higher channel thickness. It may be due to the improvement in the gate electrostatics on the channel with the decrease in tsi.

Fig. 4. Variations of threshold voltage with gate length (L) for different values of silicon film thickness tsi.

Figure 5 shows the variations of the threshold voltage with channel lengths for different values of gate-oxide thickness (tox). It can be observed from the figure that the device with lower tox shows better suppression of SCEs.

Fig. 5. Variations of threshold voltage with gate length (L) for different values of gate-oxide thickness (tox).

Figure 6 shows the variations of threshold voltage with channel length (L) for different values of stack layer thickness (tk). The threshold voltage roll-off is higher for large value of tk. It may further be noted that if we compare Fig. 5 with Fig. 6, then it may be observed that the effect of tk oversees the effect of tox on the threshold voltage variation. Even the threshold voltage roll-off is improved in the case of larger value of tk therefore in the case of TMGS-DG MOSFET gate stack oxide acts as a controlling gate oxide.

Fig. 6. Variations of threshold voltage with the gate length (L) for different values of stack layer thickness (tk).

Figure 7 shows the effects of change of high-k stack material on the threshold voltage of the device with different gate lengths. Three different dielectric materials are considered i.e., HfO2, Al2O3, and SiO2 with largest value of dielectric constant for HfO2 and least value for SiO2. It can be observed from the figure that for the material with high dielectric constant threshold voltage, the roll-off is considerably reduced. It may be due to the fact that increasing dielectric constant of the stack material can lead to a decrease in the effective gate oxide thickness teff (Eq. (1)), which may result in an efficient suppression of SCEs. This implies that the improved gate controllability can be achieved by using material with high dielectric constant to establish the suitability of TMGS-DG MOSFET over conventional DG MOSFET.

Fig. 7. Variations of threshold voltage with gate length (L) for different stack materials HfO2, Al2O3, and SiO2.
5. Conclusions

In the present paper, channel potential and threshold voltage of DG-MOSFET is modeled analytically. A detailed study to see the effects of gate stack on the threshold voltage of TMGS DG-MOSFET and other important parameters is performed in detail. The model results are found to be well matched with the ATLAS numerical simulation results for different values of drain bias, silicon film thickness, gate-oxide thickness, stack layer thickness and stack materials.

Reference
1Wong H2012Nano-CMOS Gate Dielectric EngineeringNew YorkCRC Press161–6
2Maiti CKMaiti TK2013Strain-Engineered MOSFETsNew YorkCRC Press171–7
3Ishimaru V 2008 Solid-State Electron 52 1266
4International Technology Roadmap for semiconductors 2007
5Agrawal B1994Ph. D. Dissertation, Rensselare polytech. Inst. Troy, New York
6Chaudhry AKumar M J2004IEEE Transactions499
7Chenming Hu1996International Conference on Electron Devices Meeting (IEDM 96)319322319–22
8Goel KSaxena MGupta MGupta RS 2006 IEEE Trans. Electron Dev 53 1623
9Razavi POrouji A A D2008Electronics Conference 11th International Biennial Baltic836
10Chaves F AJimenez DRuiz F JGodoy ASuñé J 2012 IEEE Trans. Electron Dev 59 2589
11Abdi M ADjeffal FMeguellati MArar D200916th IEEE International Conference on Electronics, Circuits, and Systems ProceedingsDecember 13, 2009892895892–5
12Razavi POrouji A A2008International Conference on Advances in Electronics and Micro-electronics(ENICS’08)111411–4
13Pradhan K PMohapatra S KAgarwal P KSahu P KBehera D KMishra J2013Microelectronics and Solid State Electronics21
14Alamgir M IGulib A UAhmed K MFarzana E2012International Journal of Technology Enhancements and Emerging Engineering Research196
15Li JLiu H XLi BGao LYuan B 2010 Chin. Phys 19 107301
16Solomon P Met al. 2003 IEEE Circuits and Devices Magazine 19 48
17Dao-Ming KQi LJun-Ning CShan GLei L20068th International Conference on Solid State on Integrated Cirvuit Technology (ICSICT 06)October 23–26Shanghai, China
18Young K K 1989 IEEE Trans. Electron Dev 36 399
19Oh S HMonroe DHergenrother J M 2000 IEEE Electron Dev. Lett 21 445
20Diagne BPrégaldiny FLallement CSallese J MKrummenacher F 2008 Solid-State Electron 52 99
21Chen QAgrawal BMeindl J D 2002 IEEE Trans. Electron Dev 49 1086
22Tsormpatzoglou ADimitriadis C AClerc RPananakakis GGhibaudo G 2008 IEEE Trans. Electron Dev 55 2512
23El Hamid H AIñìguez BGuitart J R 2007 IEEE Trans. Anal. Electron Dev 54 572
24Lundstrom M 1997 IEEE Electron Dev. Lett 18 361
25Jankovic N DArmstrong G A 2004 Microelectron. 35 647
26Dubey STiwari P KJit S 2010 J. Appl. Phys 108 034518